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Abstract. In this paper, we study direct and indirect Galerkin method for 
solving linear Integral- Algebraic Equations of index 1. Convergence of indirect 
method is also analyzed. 



1. Introduction 

Recently numerical solution for Integral-Algebraic Equation (IAE) has been 
studied by many researchers [1]. In physical models with constrained time or space 
variables, Algebraic Equations including Integral or Differential Equations arise. 
Therefore, the investigation of mixed types of equations is important. This can 
be seen by observing the number of papers that are about Differential-Algebraic 
Equations, Integro Differential Equations, and so on. One of these systems which 
mixes Volterra Integral Equations of first and second type is the system of Integral- 
Algebraic Equations. However, we dont find any application for this systems except 
[6, 7, 8, 9], but it seems that they be very important in the future. In this paper we 
use the standard Galerkin method, for the system of Integral Algebraic Equations 
of index 1. We implement the Galerkin method by considering the structure of IAE 
of index 1 in two ways, direct and indirect and we analyze indirect method here. 
We Consider the system 

(1.1) A{t)y{t) + [ k(t, s)y( S )ds = f(t), t E / := [0, T], 



Jo 

where A E C(I,R rxr ), f E C(I,W) and k E C(B, M rxr ) with© := {(t,s) : < s < 
t < T}. If A(t) is nonsingular for all i G / then multiplying (1.1) by Ar 1 changes it 
to a system of the second kind Volterra Integral Equations that it's theoretical and 
numerical analysis have been almost completely investigated [3, 4]. Otherwise, if 
A(t) is a singular matrix with constant rank for all t E I, then, the system (1.1) is 
an IAE. The classification of this system done by introducing index notation. The 
differential index of differential-algebraic equation can be extended to this system. 
The following extension is done by Gear in [5]. 

Definition 1.1. We say the system (1.1) has index m if, m is the minimum possible 
number of differentiating (1.1) to obtain a system of second kind Volterra Integral 
equation. 

For example the following system 

(1.2) x{t)=h{t)+( k 11 (t,s)x( S ) + k 12 {t, S )y( S )ds, 

Jo 
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(1.3) Q = f 2 (t)+ k 2L (t,s)x(s)+k 22 (t,s)y{s)ds. 

Jo 

with the functions fi, i,j = 1,2 arc sufficiently smooth and /2(C)) = 0, 

|fc 22 (M)l > k Q > for all t £ [0,T]; (1.2) and (1.3) is the system of Integral- 
Algebraic Equations of index 1. For this system the Piecewise Polynomial Collo- 
cation method is investigated by Kauthcn in [2]. In this paper we investigate the 
Galcrkin method for this system. 

2. PRELIMINARY 

The Galcrkin method is a projection type method [3, 4] that we apply to IAE 
of index 1. Let K, : X — > X be an operator, where X is a Hilbcrt space with 
orthonormal base and X„ be an n-dimensional subspacc of X generated by 

{0i}" =1 . Then the system (1.2) and (1.3) changes to 

(2.1) Kx = f. 

For using the Galerkin method's we define the projection P n : X — > X„ by P n x = 
Y^i=i( x > Then it is clear that ||P ra x — x|| — > as n — > 00. We search x„ £ X„, 
x„ = xcf) with x = (xi , . . . , x n ) and <f) = (<pi , . . . , <ft n ) T , such that 

(2.2) P„/Cx„ = P n f. 

The corresponding operator of the system (1.1) is defined by 
(2.3) 

K{\x,y)) = (x- j k 11 (t,s)x(s) + k 12 {t,s)y(s)ds, - J k 21 (t,s)x{s) + k 22 (t,s)y(s)ds 

and the assumption that / = (/i,/2) T , X C (£ 2 [Q,T]) and the functions fi, fc u , 
i,j = 1,2 are sufficiently smooth and /2(C)) = 0, \k 22 (t,t)\ > k a > for all t £ 
[0, T], guarantee us that K,x = f has a continuous solution as it follows. We first 
differentiate equation (1.3) w. r. to t and we have 

=f' 2 (t) + k 21 (t,t)x(t) + k 22 (t,t)y(t) 

* dk 21 (t,s) . . dk 22 {t,s) . _ 

Because \k 22 (t,t)\ > fc > we obtain 

k 22 {t,t) k 2A (t,t) 

1 dk 21 {t, s ) . , 1 dk 22 {t, s ) , , , 
-< s ) - u22n * — m — y( s ) ds ' 



lo k 22 {t,t) dt v ' k 22 {t,t) dt 

as a Volterra integral equations of the second kind for y(t). Now it is easy to show 
that this system has a unique solution (see [4]). 

3. Implementation and Convergence 

We can implement the Galcrkin method for two systems (1.2), (1.3) and (1.2), (2. 4). 
For the former the Galerkin method is called direct and for the further, is called 
indirect method. For the second system, the convergence and order of convergence 
can be obtained by similar arguments as discussed [6] and [3]. In [3] the order of 
convergence has been explained for Fredholm integral equation. Since, by using 
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suitable kernel, each Volterra integral equation can be changed to a Fredholm in- 
tegral equation, Galerkin method for both are the same. Hence the convergence 
order of the approximate solution is equal to the convergence order of the best 
approximation in space X n . This is confirmed by numerical examples. Let {1^}°^ 
be Legendre polynomials as orthonormal base for €?[— 1, 1]. For £ 2 [0,T] we use 

shifted Legendre polynomials {FtJ^i, V Tx {s) = ^V.i 2 ^^), s e [0, T] and for 

-2r n Tl^2 



(£ 2 [0,T]) 2 , we use {(V Ti , 0), (0, V*)}^ with ({x,y), (z,t)) = (x,z) 2 + (y,t) 2 and 
{x,y)2 = J Q xydt. Also we suppose P 
equations (1.2), (1.3) and (2.4) we have. 



x,y)2 = Jq xydt. Also we suppose P n ([x,y]) = (Yn=i x i V Ti'YH=iyi v Ti) ■ For 



(3.1) 



" n „t 

i=i i=i J ° 

n i-t 

+ X> / k 12 (t,s)V Tl (s)ds, 
<=i Jo 



t 

22 



(3.2) = Mt)+J2xi k 21 (t, S )V T Mds + Y,y* k 



(t,s)V Ti (a)ds, 



^ "g^/o fe 22 (M) ft VTMdS 

-j>/ p 2 (m) ^ ^ (s)ds - 

Now we can use (3.1), (3. 2) for the first system and (3.1), (3. 3) for the second system 
introduced above. It is important for analysis to note that formula (3.3) is derivative 
of (3.2). Now if we multiply equations (3.1)-(3.3) to Vr- (t) and integrate on the 
[0, T], we have obtained 



r-l " />T r-t 

: / fi(t)V Tj (t)dt + J2 x i / / k ll {t, S )VTAs)V T] {t)d S dt 
Jo i=1 Jo Jo 

n rT rt 

+ Y> / / k 12 (t,s)V Ti (s)V Tj (t)dsdt, 

Jo Jo 



■T, 



(3.4) 



0=1 h(t)V T] (t)dt + Vx,/ f fc 21 (i, s)V Ti ( S )V T] (t)dsdt 
Jo i=1 Jo Jo 

+ Y,yi / k 22 {t,s)V Tz {s)V T] (t)dsdt, 
„•_ 1 Jo Jo 



(3.5) 
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(3-6) -E*/ ^'feW^W* 

The direct and indirect Galcrkin method for system (1.2), (1.3) are solving the linear 
systems (3. 4), (3. 5) and (3. 4), (3. 6) respectively. 

Lemma 3.1. [4] Let V and W be normed spaces, with W complete. Let K G 
C(V,W), let {K n } be a sequence of compact operators in C(V,W), and assume 
K n — > K in C{V, W). Then K is compact. 

Lemma 3.2. [4] Let V be a Banach space, and let {P n } be a family of bounded 
projections on V with 

P n u — > u as n — > oo, u G V 

If K : V — > V is compact, then 

\\K-P n K\\->0 as n^oo. 

Theorem 3.3. [4] Assume K : V — » V is bounded, with V a Banach space; and 
assume A — K : V — > V is one to one and surjective. Further assume 

\\K - P n K\\ -> as n -> oo. 

Then for all sufficiently large n, say, n > N, the operator (A — A') -1 exists as a 
bounded operator from V to V. Moreover, it is uniformly bounded: 

sup || (A - P n Ky 1 \\ < oo. 

n>N 

For the solutions u n (n sufficiently large) of P n (X — K)u n = P n F and the solutions 
u of (A — K)u = F, we have 

u — u n = A(A - P„A') _1 (u - P n u) 

and the two-sided error estimate 

l A l IL. o ..II ^ II.. .. II ^ IM \\t\ D 



— ||u-P n u|| < ||«-u„|| < |A| \\{X-P n K)-'\\ ||u-P n u|| 



\\X-PnK\ 

We have following Theorem for indirect method. 

Theorem 3.4. Let for the system (1.2) and (1.3) the following conditions are 
satisfied: 

(1) ft G C 1 (0 > T) > KU[t,8] G C^B), for alli,j = {1,2}, 

(2) k 22 (t,t) > k > 0, 

where D := {(7;, s) : < s < t < T}. Then the convergence order of the approximate 
solution is equal to the convergence order of the best approximation in the space 

Xu = {VtA°Zi- 
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Proof. The space £ 2 [0,T] with C 2 Norm is Banach space, hence the Cartesian 
product on this space, i.e., V = (£ 2 [0,T]) 2 , with introduced norm generate also 
Banach space. Now we change the system (1.2) and (2.4) to the following system 

(I-K)x = F, F=[F 1 ,F 2 ] T , x= [x u x 2 ] T , 



(3.7) 
with 



K (t, s) 



F 1 (t) = f 1 (t), F 2 {t) 
fc u (M), s<t, 



f 2 (t) k 2 \t,t) 



k 22 (t,t) k 22 (t,t) 



AW. 



0. 



s > t, 



K i2 (t,s) 



k 12 (t,s), s<t, 
s > t, 



K 21 (t,s) 



0. 



k 2i (t,t) K V> s ) 



0. 

dk 21 (t,s) 



dt 



K 22 (t,s) = 



0. 



k 21 (t,t) j 12( , \ 



k 22 (t,t) 



dk 22 {t,s) 



dt 



and with the integral operator 

r K n (t,s) K 12 (t,s) 
K 21 (t,s) K 22 (t,s) 



Kx = 
By introducing 

K„x := 

with continuous kernels 
K?(t,s) = 

K 12 {t,s)-- 



K?(t,a) K} 2 {t lS ) 
Kl\t,s) Kl 2 {t lS ) 



xi{s) 
x 2 {s) 



xi(s) 
x 2 (s) 



s<t, 
s>t, 

s<t, 

s > t, 



ds. 



ds, 



fc n (M), 



s < t, 



-n(fc u (M))((* +£)-*), se[t,t+i], 

o, s >t + b 

k 12 {t,s), s<t 

-n(ki 2 (t,t))((t + ±)-s), se[t,t+±}, 

0, «>t+i 



_n (fc22(* t 't) fcll ( i )i) T fei!2( tit ) 



1 dk* L {t,s) 
k 22 {t,t) dt ' 

1 ^^((t +£)"-). 



0, 



~fc^(t,t) 



1 dk 22 {t,s) 
k 22 (t,t) dt 



s <t, 
s < i, 



0, 



p^i. (M) + F ^y^F M )((t + i)- S ), -€[t,t+i], 

s >t+i. 



which is compact and using lemma 3.1, imply K is compact. 

Applying the Galerkin method to the Frcdholm system (3.7) leads to the linear 
system equivalent with the system (3. 4), (3. 6). So it is enough to show that for this 
system the convergence order of the approximate solution by the Galerkin method 
is equal to the convergence order of the best approximation. Since, the system 
(3.7) also is a volterra system, using the same line of [4], (/ — K) is invcrtible so is 
surjective. The Lemma 3.2 provide the remainder condition. Finally all conditions 
of the theorem 3.3 are satisfied and this completes the proof. □ 
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Table 1. Numerical result obtained from (3. 4), (3. 5) 



n 


2 


4 


6 


8 


10 


\\x n - X\\ 


4.0e - 002 


3.0e - 004 


7.4e - 007 


9.4e - 010 


7.6e - 013 


Wvn - v\\ 


1.6e - 001 


2.2e-003 


8.1e-006 


lAe - 008 


1.4e-011 



Table 2. Numerical result obtained from (3. 4), (3. 6) 



n 


2 


4 


6 


8 


10 


\\x„ - x\\ 


3.8e-002 


2.6e - 004 


6.7e - 007 


8.8e-010 


6.8e-013 


hn - y\\ 


7.3e - 002 


5.1e - 004 


1.3e - 006 


1.7e - 009 


1.3e - 012 



Table 3. Logarithm of the error by best approximation using Lc- 
gcndrc polynomial) 



n 


2 


4 


6 


8 


10 


||P„sin(t) -sin(t)|| 


5.1e - 002 


3.3e - 004 


8.1e-007 


l.Oe - 009 


8.0e-013 


||P„cos(t)-cos(f)|| 


7.8e - 002 


5.4e - 004 


1.4e - 006 


1.7e - 009 


1.3e - 012 




-14 1 ' • • • • • • 1 -12 1 1 1 1 1 1 1 1 1 

2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 10 



FIGURE 1 . plots of convergence order by direct Galerkin method 

4. Numerical results 

In this section we illustrate the methods by numerical examples. 
Example 1. 

We consider (2.3) with k 11 = s + 1, k 12 = s 2 + t 2 , k 21 = s - 1 2 , fc 31 = s + 1 + 1, /i = 
-t-2sin(t)t 2 + 2sin(t) and / 2 = t 2 ~2sin(t)+cos(t)t-cos(t)t 2 + l-cos(t)-2sin(t)t 
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FIGURE 3. plots of convergence order of best the approximation 



where the exact solutions are x = sin(t) and y = cos(t). Numerical results obtained 
by two methods are in table 1 and 2 which show rapid convergence. In these tables 
ll/H = max T6 [ 0jT ] \f(t)\. In the following figures, we also illustrate log 10 (||x„ — x\\) 
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and log 10 (||j/ n — y\\) . For the second system, we expect that this figure to be similar 
to the figure of log 10 (|jP n x — x||) and log 10 (|jP n ?/ — y\\). This numerical experiment 
shows that the second system is more accurate than the first. But we need more 
computations in comparison with the first. 

Remark 4.1. In equations (3. 4), (3. 5) and (3.6), we find two type of integrals that 
should be computed, single and double integrals. The integral of type L f(t)dt 
can be computed by Gaussian Integration Methods (see [11]) of the form 

and by changing of variable we have 

(4.1) / f{t)dt~ J2wif(Txi/2 + T/2) 

where the roots of Vr, n +i)i an d w i, — i w n are computed by the system 

expressed in (3.6.13) of [11]. An approximate values of are computed as follows 
So So f^ti s )dsdt in this calculation is as it follows 

(4.2) / / f(t, s )d S dt = Y,Y, w wtt^r + ^ 2 t 2 +- L i- L ) 

JO Jo i=li=l III I 

Also we notice that the term n is not fixed for different P m , indeed in our calculation 
ti = tyl that means for improving accuracy of the solution we improve the accurate 
approximate solutions of integrals. 
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